DETERMINING THE ORDER OF 5*7(3) 
DECONFINING PHASE TRANSITION 

Zheng-Kun Zhu 

Department of Physics and Atmospheric Science, Drexel University^ 

Philadelphia, Pennsyvalnia 19104-9984 
Center for Theoretical Physics, Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 
and 

Da Hsuan Feng 

Department of Physics and Atmospheric Science, Drexel University 
Philadelphia, Pennsyvalnia 19104-9984 
Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37996 

Abstract 

An effective spin modei for the finite temperature non-abelian lat- 
tice gauge theory is derived. The outcome is the surprising result 
that only nearest neighbor coupling survives, thus confirming the well 
known numerical results that the deconfining phase transition of the 
(3+l)-dimensional SU(3) pure gauge theory is first order. 
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A major success of lattice gauge theory (LGT) is the numerical demon- 
stration of a transition from the low temperature color confining to the high 
temperature non-confining phase]]]]. Determining the order of the transition 
is a subtle matter. For the (3+1) dimensional 577(3) finite temperature LGT, 
it was confirmed numerically on the N T = 4 lattice (N T is the lattice size in 
the temporal direction) that the transition is first order 0, |3|, |J. More recent 
calculations on the N T = 6@] also suggest that this is the case. Still, the 
definitive statement about the order of phase transition remains a challenge 
for the next generation of computers. 

Apart from the full numerical determination the order of the finite tem- 
perature phase transition, there exists in the literature universality argu- 
ments, first suggested by Svetitsky and Yaffe in 19820, which relate the 
field theory to a simpler three dimensional effective spin model. The deriva- 
tion of the effective spin model is highly complex. Nevertheless, by conjec- 
turing a short range coupling, it was shown that the transition is also first 
order ||, [7]]. Although the conjecture may not be beyond reproach, these 
important results certainly lend confidence to the conclusion that the tran- 
sition could indeed be first order. We therefore feel that there remain the 
urgent task of deriving an effective spin model. The purpose of this letter is 
to carry out the derivation. 

Instead of directly studying the (3+1) dimensional SU(3) finite temper- 
ature LGT, we will begin by studying a generic (3 + l)-dimensional SU (N) 
LGT. The finite temperature behavior can be described by a partition func- 
tion defined on a hypercubic lattice of size Ng x N T , 

Q = Ee 5 , (1) 

where S is the Wilson action 

5 = &EE ^ReTr{UltK +ll UCK + ), (2) 

(3e is the coupling constant, n and /i, v represent the space-time coordinate 
and directions respectively. Finite temperature is introduced by imposing a 
periodic boundary condition in the time direction, with period N T . Accord- 
ingly, the temperature is jr-, where a is the lattice spacing. 

It is well known that the nature of the deconfining phase transition can 
straightforwardly be studied by constructing an effective theory (of the Wil- 
son line) from the partition function in eq.([L]) with all the link variables 
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integrated out except those for the Wilson lines. Then the effective action 
SYff has the following form||: 

S%f = S% f ({^TrW n }). (3) 

Unfortunately, this process could only be carried out in the strong coupling 
limit. 

There are two main obstacles to obtain S^t. The first is that in eq.(^) 
there is no explicit Wilson line variable in S. The second is that it is very 
complicated to carry out the integration. We will now discuss the removal 
of these difficulties. 

It turns out that by taking a thermal gauge choice0 the first difficulty 
can easily be removed. 

XJl^ = l(l<n T <N T -l) 

The link U^ Nt remains unchanged, and the trace of U^ Nt becomes a Wilson 
line, relabelled here as W n . By inserting the above gauge choice into eq.(||), 
we can rewrite the action S of eq.(fj) as a sum of S 9 and S T , where 

S 9 = As E E ^Tr(Ul n Ui +i>n U^ +jjn X% T ) 

n,i<j n T =l 

N T -1 1 

E JrReTr(Ul n U£,n T+ i) (4) 

n,i n T =l 

and 

^ = f3EY.l ReTr KNW Q+t u£ 1 W+). (5) 

n,i 

Straightforward removal of the second difficulty is arduous. To this end, 
we will instead effectively decouple the partition function into two indepen- 
dent sub-partition functions; one describes the Wilson line field and the other 
the space-like link field. We will then derive the effective spin model by means 
of a variational principle. In the following, we will show how to decouple the 
partition function. 
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First we will formally write the effective action as a sum of two parts: one 
is an effective spin model, the other has only space-like link variables ||f0||. 



Se f f = SY ff ({W n })+SV ff ({Ui}). (6) 



where Sfff and S^t are 



Q w , ^{u^} ex P(S 9 + S T ) 



'eff 



2^{W} eX PWeff 



We can now see that the partition function of the effective action given by 
eq.(f|) is the same as that of the action given by eqs.(||) and 

» S " ( 9 ) 

u,w u,w 

In this way at least the decoupling of the fields W and U can formally be 
achieved. We can now see that W and U fields can be described exactly by 
SYff and Sgft respectively. It should be noted that S e /f cannot describe the 
correlation between W and U. However, such a correlation is irrelavant for 
the physics discussed in this letter. 

Of course, S^j can in principle be obtained by intergating U in eq.(0). 
However as we have discussed earlier, it is highly complex to obtain S£ft by 
this approach. Instead we have resorted to the variational principle to derive 
SYff- Since S^,* does not explicitly connect with the deconfining transition 
in our scheme, we will not comment on it further. Indeed, there is no loss of 
generality if we simply assume that it is already known. 

The action S e ff resembles S except for the coupling between W and 
U fields. Therefore it is quite reliable to determine S^f by a variational 
method|TT| where S e ff is a trial action. To this end, we will first determine 



its form. We notice that there is a local gauge invariance for W in S of 

eqs.(f|) and Hence, in order to maintain this invariance, SYff must only 

depend on Tr(W m ) (where m is an integer) Jl2| . Then we can describe the 
<w 

eff 



action S3 f as follows: 



Sfff = PEaY,M^ T TrW n+ ~TrW+) + S?({±Tr(W?)}). (10) 



' N N n/ N 
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where a is a variational parameter. From eq.([|), we note that the first term 
of eq.( |IO|) is the dominant part of the effective spin model and the residue. 
What remains is to determine a and the form of . 

Ignoring the residue part in eq.(|10D for the moment, we write the action 

as: 

= a(3 E J2Re(^TrW n+ ~TrW+), (11) 
then the trial action becomes 

So = S + S e ff. (12) 
To determine the value of a, we calculate the partition function as follows, 
Q = Q < exp{S 9 + S T - S } > , (13) 
where < • • • > represents the average in action So and Q is 

Qo = £e s °. (14) 
Using Jensen's inequality (< expX >> exp < X >), we can obtain 

F = lnQ> lnQ + < S 9 + S T - S > (15) 
and a can then be determined by maximizing F var with respect to a 

F mr = lnQ + <S 9 + S T -S >o- (16) 
To accomplish this, we will compute < S 9 + S T — So >o ( = K)\ 

^ = 3(p E Re(<^TrW+^TrW n+l > W o<^Tr(Ul N Ul + tl )>^ 

-a(3 E Re(< ^TrW+ ^TrW n+i > wo )) 
<S 9 - S%f >y 

where < • ■ ■ >u and < ■ • ■ >wo represent the average in S^tf and S^ respec- 
tively. We then notice that S^ is real and invariant under the transformation 

W -> W + , (18) 
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hence < -^TrW^ j^TrW n+ i >wo is also real. By maximizing F var with re- 
spect to a , we obtain 

a = Re{< ^Tr(Ul N U^) >u) (19) 
We are now ready to determine the form of Let's first write as 

= ^ M (20) 

where £ is a variational parameter and S^ M the remaining part of the effec- 
tive action which is independent of £. 

We notice that the action S of eqs.(f|) and @ is invariant if all the link 
variables U transform as U — > U + (here U means either U or W). To maintain 
this invariance in the effective spin model, S^f must be invariant under the 
transformation W — ► W + . 

Since SYff is real and invariant under W — ► < -^TrW+^TrWn+j >^ 
(where < ■ ■ • >iy represents the average in action S^j,) is still real. With 
this, we will again compute the partition function Q e ff in the action S e ff 
and K as: 

K = <S 9 -S^ ff >u-<S^> W , (21) 
By maximizing F var = lnQ e ff + K with respect to £, we obtain 

gg^ a < S™ > w _ n , 99 s 

"af — " ( } 



From eqs.(Efl) and (|22|), we obtain 



and 



< >W = 0. (23) 



a ^ cWM 

^ ' ; > ^ =< (^ m ) 2 >w - < sr M >^ (24) 



Clearly according to eq. (j24"|) , 9<Sr d ^ >w does not vanish unless S^ M is a 
constant, a situation which is not relavant in our discussion. From eq.(|23|), 
we see that £ must vanish. Our derivation is general since S^ M includes any 
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kind of coupling terms. Hence surprisingly, we have obtained an effective 
spin model with only nearest neighbor coupling terms. 

It is important to understand why by using the variational principle, 
the effective spin model appears simple. Actually, it is very complex in the 
sense that to compute a according to eq . (|19D is complicated. The reason is 
as follows: In order to obtain a, we have to know the form of Sgft which 
can be obtained by integrating over W in eq.(§). With this knowledge, a 
can then be calculated according to eq.(^). Both procedures are nontrivial. 
In fact, the difficulty of deriving the effective spin model has actually been 
transferred from the integration of U in eq.(p]) to the integration of W in 
eq.@ and to calculate a from eq . (|19D . Therefore the simplicity of the form 
of the effective spin model is not obtained without a price. What we have 
done is to "organize" the theory such that the discussion of the order of 
deconfining phase transition can be vividly studied. 

Having the effective spin model, we can now focus on studying the de- 
confining phase transition. To this end, we note that there were numerous 
investigation made in this direction already. These studies enable us to dis- 
cuss the order of the deconfining phase transition. 

For the SU(2) gauge theory, it is well known that the plaquette is con- 
tinuous with the change of 0e- Therefore, according to eq.fjji]), a is also a 
continuous function of (3e- Hence our study directly shows that the phase 
transition of the Z{2) spin model possesses the same universality property 
as the SU{2) gauge theory. This agrees with the recently reached conclusion 



by the Monte-Carlo real space renormalization study[13|, 14 . 

For the SU(3) gauge theory, it was shown via the Monte Carlo study of 
the effective spin model with a nearest neighbor coupling assumption that 
the transition is first order [l], |7|. Results reported here confirm the nearest 
neighbor assumption and are consistent with the recent Monte-Carlo results 
about the transition order 0, 0. 

Finally, we mention that our derivation appears to provide an analytical 
approach to study the thermodynamical behavior of LGT|15| . 
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